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Abstract 

Hadron transverse momentum spectra are proposed as a means to probe the under- 
lying partonic dynamics in deep inelastic scattering. The BFKL evolution equation, 
postulated for small Bjorken-x, leads to an enhanced parton emission over the con- 
ventional DGLAP ansatz, and can thus be tested. 



1 Introduction 



The successful description of the nucleon structure function data by perturbative QCD, 
cast into the DGLAP (Dokshitzer-Gribov-Lipatov-Altarelli-Parisi) parton evolution equa- 
tions [I] constitutes one of the major successes of QCD. At small enough Bjorken- 
x however, these equations are expected to break down. An alternative ansatz for the 
small x regime is the BFKL (Balitsky-Fadin-Kuraev-Lipatov) equation 0. At lowest 
order the BFKL and DGLAP equations resum the leading logarithmic (a s lnl/x)™ or 
(a s ln(Q 2 /Qq)) u contributions respectively, with Q 2 being the virtuality of the exchanged 
photon. The leading log DGLAP ansatz corresponds to a strong ordering (Qq -C k T \ <C 
...kj- 2 -C ...Q 2 ) of the transverse momenta kj- (w.r.t. the proton beam) in the parton 
cascade (Fig. [1|) while in the BFKL ansatz they rather follow a kind of random walk 
[k T 2 ~ kx 2 +l ) ||. It is an open theoretical question to what extent BFKL type con- 
tributions play a role in the small x (~ 10~ 4 ) regime now accessible in deep inelastic 
scattering (DIS) at HERA. Measurements on the hadronic final state emerging from the 




Figure 1: Parton evolution in the ladder approximation. The longitudinal fractional momenta Xi and 
transverse momenta kri of subsequently emitted partons are indicated. 

cascade could be sensitive to the new type of evolution. For example, without the re- 
striction of strong ^ ordering, more transverse energy Et is expected from BFKL than 
from DGLAP parton radiation in a region between the current region and the proton 
remnant || . Though the HERA Et flow data || can be interpreted consistently with the 
BFKL mechanism, it was not possible to disentangle the perturbative parton radiation 
from non-perturbative hadronization effects || [jj] . 

In this paper it will be demonstrated that single particle transverse momentum (pt ) 
spectra represent a more direct measure of the partonic activity than the Et flow mea- 
surements. Observables are then constructed which allow to discriminate the fey ordered 
from the unordered parton shower scenario. 

Predictions for the cases of the ordered resp. unordered cascades are extracted from 
Monte Carlo models, which incorporate the QCD evolution in different approximations 
and utilize phenomenological models for the non-perturbative hadronization phase. The 
MEPS model (Matrix Element plus Parton Shower) ||, incorporates the QCD matrix 
elements up to first order, with additional soft emissions generated by adding leading log 
parton showers. In the colour dipole model (CDM) |], [K| radiation stems from colour 
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dipoles formed by the colour charges. Both programs use the Lund string model [O for 



hadronization. The Herwig model JT2[ is also based upon leading log parton showers, with 
additional matrix element corrections 
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It uses a cluster fragmentation scheme 
The CDM description of gluon emission is similar to that of the BFKL evolution, because 
the gluons emitted by the dipoles do not obey strong hr ordering [15]. In MEPS and 



Herwig the partons are strongly ordered in hp , because they are based upon leading log 
DGLAP parton showers. The latest versions of the models (Lepto 6.4 for MEPS, Ariadne 
4.08 for CDM and Herwig 5.8) are used with the parton density parametrization MRSH 
||16|| . They provide a satisfactory overall description of current DIS final state data [17| , 
in particular of the Et flows Q. 



2 The Method 

In this section the sensitivity of single particle pr spectra to the parton activity in the 
ladder is demonstrated. Generated events are selected from the kinematic plane of x and 
Q 2 according to the binning chosen by HI Events from two bins, one at "low x " 
({x )=0.00037) and one at "highx " ({x )=0.0023), with (Q 2 ) « 14 GeV 2 approximately 
constant, are compared. In Fig. ^| a the Et flow in the hadronic centre of mass system 
CMS is shown Q as a function of pseudorapidity r\ [r\ = — In tan 9/2, where the angle 9 
is measured w.r.t the virtual photon direction) for events with small x . As expected, 
the partons produced from unordered emission (CDM) give more Et in the central rj 
region rj « than the ones emitted from the ordered cascade (MEPS, Herwig). However, 
the observable particles emerging after hadronization give rise to very similar Et flows, 
unresolvable with current data [f|, [17| . While hadronization adds relatively little Et to 
the partonic Et for CDM, most of the Et is generated by hadronization in the cases 
of MEPS and Herwig. To answer the question whether the Et observed in the data is 
generated predominantly by parton radiation or by hadronization, inclusive pt spectra 
are considered. Hadronization should produce typical spectra which are limited in pt , 
while parton radiation should manifest itself in a hard tail of the pt distribution. That 
tail is due to occasional hard parton radiation, from which hard particles can emerge. 
The production of such hard particles from hadronization would be suppressed. 

To test this idea particles from a "central" r] interval < rj < 2 are examined. The 
lower limit is given by the approximate acceptance of the HERA detectors, and the upper 
limit restricts the interval to the region where the partonic differences in Et are largest, 
excluding the "current" fragmentation region. Events are compared which have similar 
hadronic Et in that interval [E^ between 1 and 2 GeV/unit rapidity), but different 
amounts of partonic Et , E T ar ■ Events with E T &r < 0.2 GeV/unit rap. are called 
hadronization dominated, and events with E T W / > 0.5 are called parton dominated. 
The correlation between E T ai and Et is shown in Fig. b. For the CDM two classes 
of events can be identified. For one class E T ad is well correlated with E^ , for the other 

1 In MEPS the new concept of soft colour interactions (l]] had to be introduced to reach the level 
of Et seen in the data [Q, [TtJ . Intriguingly, this mechanism also produces rapidity gap events jl8) at a 
rate comparable to observation 0, roughly 10%. Rapidity gap events are also produced by the cluster 
fragmentation in Herwig. In this paper rapidity gap events are excluded. 

2 All distributions shown arc normalized to the number of events N which enter the distribution. 
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Figure 2: a) Transverse energy flows vs. pseudorapidity rj generated by the models CDM, MEPS 
and Herwig for partons and for hadrons at "low x ". The proton direction is to the left, b) Correlation 
between partonic and hadronic Et produced in the central pseudorapidity bin < r\ < 2 for CDM, 
MEPS and Herwig. 



Ej? r is small, regardless of E T ad . For the other models, most events fall into the latter 
class, while the correlation between E^ r and is much less pronounced for the rest 
of the events. The parton dominated events indeed exhibit a harder p? spectrum than 
the hadronization dominated events (s. Fig. |3] a), regardless of the underlying parton 
dynamics or the applied hadronization model. Therefore px spectra provide a useful 
method to study the underlying parton dynamics in DIS. 



3 Predictions 

In this section observables are constructed that should allow to distinguish between the 
two scenarios of ordered resp. unordered parton evolution, or in general be sensitive to 
the parton radiation generated in the evolution. In this study the CDM is taken as a 
model to represent the unordered parton cascade, and the MEPS and Herwig models 
represent the ordered cascade. In Fig. |3] b the inclusive px spectra of charged particles 
from the "central" rj bin are shown for large and for small x . At large x , all models 
predict similar px spectra. At small x however the tail of the distribution (p? larger than 
~ 1.5 GeV ) is harder for the unordered model (CDM) than for the others, as expected 
given the larger parton activity. Less visible due to the logarithmic scale is a difference 
in the average charged multiplicity in the central rj region between the two scenarios of 
about 20%. In MEPS and Herwig more soft particles are produced in the hadronization 
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Figure 3: a) Transv. momentum (j>t ) spectra of charged particles from the central pseudorapidity 
bin < rj < 2 for events which are either hadronization dominated or parton dominated, b) Ch. particle 
Pt spectra from < rj < 2 for events with "high x " and with "low x ". c) Event distribution in 
the quantity E^ d , determined from < ?y < 2, at "low x ". d) Distribution of the maximal transv. 
momentum p™ ax of ch. particles from < rj < 2 for events with E^ d > 3 GeV /unit rap. at "low x ". 



phase to generate the E T seen in the data. 

One also notices from Fig. |3| b that from the unordered cascade a hardening of the 
spectrum is predicted towards small x , and a softening otherwise. This behaviour can 
be traced to the fact that while all models predict an increase of E^ d towards small x , 
as observed in the data [|], only the CDM shows that increase also on the parton level 
(see Fig. |]). The other models predict a decreasing . This behaviour of the models 
is in accord with perturbative calculations of the central Et as a function of x ||, based 
upon either the BFKL or DGLAP evolutions. As a consequence, the relative amount of 
hadronization to E^ ad decreases for CDM, but increases for MEPS and Herwig towards 
small x . The signal can be enhanced by selecting events in which large £ , ^ ad (s. Fig. || c) 
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Figure 4: Et flows vs. pseudorapidity for CDM, MEPS and Herwig on the parton and on the hadron 
level. Compared are the Et flows for "high" and "low" x . 

is observed (e.g. calorimetrically) , because in CDM that is correlated with large E T aT , 
as opposed to the other models. For such events a dramatic signal can be obtained by 
measuring the maximal pr observed in the central rj region, s. Fig. |3] d. Enhanced parton 
radiation would also be signaled in the tail of the E T ad distribution, s. Fig |3] c. 

4 Conclusions 

In order to investigate the dynamical features of parton evolution in the proton at small 
x , observables based on single particle px spectra have been constructed. It has been 
demonstrated for all models investigated that the hardness of such spectra is sensitive 
to parton radiation from the cascade. Since for small enough x it is expected that the 
DGLAP equations with strong kx ordering for parton radiation cease to be valid, and may 
possibly be substituted by the BFKL ansatz, predictions are obtained for the two scenarios 
of ordered resp. unordered cascades. They have been derived from different Monte Carlo 
models which either obey kx ordering or do not underlie such a restriction. The unordered 
scenario gives rise to a harder px spectrum in the central rapidity region of the hadronic 
CMS than the ordered one. It is further predicted that the px spectrum becomes harder 
resp. softer with decreasing x for the unordered resp. ordered scenario. The application of 
the presented method at HERA would not only allow to discriminate between the different 
QCD models, it would also offer the possibility to resolve the question of kx ordered vs. 
unordered cascade, or DGLAP- vs. BFKL- like evolution at small x . 
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